Programmable networks for quantum algorithms 
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The implementation of a quantum computer requires the realization of a large number of iV-qubit 
unitary operations which represent the possible oracles or which are part of the quantum algorithm. 
Until now there are no standard ways to uniformly generate whole classes of 7V-qubit gates. We have 
developed a method to generate arbitrary controlled phase shift operations with a single network 
of one-qubit and two-qubit operations. This kind of network can be adapted to various physical 
implementations of quantum computing and is suitable to realize the Deutsch-Jozsa algorithm as 
well as Grover's search algorithm. 



The experimental implementation of complex A^-qubit 
operations (where N > 3) and the realization of com- 
plete quantum algorithms are major challenges in the 
field of quantum computation. Any progress in this di- 
rection proves the practical feasibility of quantum com- 
putation. Further, it provides a tool for systematic study 
of the physical und technological requirements for quan- 
tum computers such as parametric constraints of a given 
implementation, decoherence times, robustness of many- 
qubit entanglement, measurement efficiency etc. 

The attempts of practical implementation focus es- 
sentially on Shor's algorithm |l|, Grover's database 
search Q, and the Deutsch-Jozsa algorithm 0, 0- 
Qubit-based experimental realizations of quantum algo- 
rithms have been achieved with liquid-state NMR tech- 
niques 0, El 0, IE an d with trapped ions 0] . In order 
to practically implement (or rather simulate) a quantum 
computer one needs to complete the following steps: pre- 
pare the initial state, apply an A-bit unitary operation 
which encodes the properties of a certain function / (the 
so-called oracle), perform another sequence of operations 
(the quantum algorithm that extracts the properties of 
/) , and finally measure the qubit register which contains 
the desired information about /. 

The universality of quantum computation implies that 
it is possible to generate arbitrary A^-qubit gates by 
using se quen ces of one-qubit and two-qubit operations 
only |l ll Il2l Il3|. Barenco et al, and later Cleve et al., 
have developed methods to design networks for A^-qubit 
controlled operations [Til Hij|. While these methods, in 
principle, are sufficient to generate any A^-qubit gate re- 
quired for the known quantum algorithms, it is not ob- 
vious how to apply them for a systematic and efficient 
practical realization of the operations. 

A feature common to all existing implementations so 
far is that the sequences of one-qubit and two-qubit op- 
erations to realize the algorithm depend on the specific 
physical system and the number of qubits. More impor- 
tantly, they even depend on the particular choice among 
the possible functions /. For example, all A^ = 3 im- 
plementations of the Deutsch-Jozsa algorithm 0, 0, 0] 
use some classification of the (2 N )\/((2 N ~ 1 )\) 2 balanced 



functions and give prescriptions how to realize the func- 
tions in each class. While for N = 3 there are only 70 
balanced functions, this approach appears hard to ex- 
tend even to A^ = 4. This situation is not satisfactory. 
Clearly, scalability is a requirement not only for quantum 
hardware |l7j . but also for quantum software. 

In practical realizations of quantum information pro- 
cessing, one has to cope also with other problems: firstly, 
it is often not possible to directly perform two-qubit op- 
erations between arbitrary pairs of bits. Secondly, the 
controlled-NOT (cnot) gate is not the genuine two-qubit 
gate for many proposed qubit systems. Different inter- 
action Hamiltonians provide different types of two-qubit 
gates in a "natural way" , i.e., gates which can be achieved 
with a single two-qubit operation 0, 0, |2£j ■ 

In order to overcome these difficulties, we have worked 
out a systematic approach to design quantum networks 
which are suitable to perform arbitrary controlled phase- 
shift operations 

U r . |x> — e-^lx) (1) 

on A^ qubits (here |x) = \xi, . . . , xn), Xj £ {0,1}, de- 
notes an element of the A^-qubit computational basis). 
The most important feature of these networks is that the 
parameters 8 X of the A^-bit operations are determined 
exclusively by the rotation angles of single-qubit oper- 
ations. Therefore, these networks may be regarded as 
programmable (in a technical sense): one can imagine 
to load network instructions into the processor with a 
punch tape where the instructions are encoded by the 
2 N — 1 rotation angles. This scheme resembles the struc- 
ture of information processing in conventional computer 
processors and may be considered as a step towards a 
von Neumann Architecture for quantum computers. 

We show that this kind of network can be extended 
recursively to an arbitrary number of qubits assuming 
only nearest-neighbor coupling. As the method utilizes 
only one type of two-qubit gate, it can be, in principle, 
adapted to a large class of hardware implementations. As 
an example we discuss the implementation of the network 
with up to four Josephson charge qubits. 
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The elementary one-qubit gate in our networks is the 
phase shift (z rotation) 
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Note that — [(j)}^ does not mix the states {|0), |1)} of the 
computational basis, but generates a phase shift whose 
sign depends on the qubit state. The fundamental two- 
qubit gate is the cnot 
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To motivate the discussion of the general method be- 
low, let us consider the simple two-bit network in Fig. ^ 
The easiest way to understand this network is to ana- 
lyze its action on the states of the computational basis 
x) = \xi,x2). We see that the network leaves the ba- 
sis states unchanged, it merely generates a phase factor 
that depends on the index x of the basis state (see the 
caption of Fig. Therefore, the action on a superposi- 
tion of basis states is that the modulus of the amplitude 
for each component remains the same while the relative 
phases change in a well-defined way. 

The total phase shift # x due to the action of the net- 



work |x) 
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x) can be written more formally as 



Here we have added a (physically irrelevant) global phase 
4>oo/2 which applies uniformly to all basis states. 



FIG. 1: Phase-shifting network for two qubits. The first two 
z rotations change only the phase of the input state \x\,x-i): 

the gate — [(jJioJ^ on ^ e Xl ^ ne s ^i^ s ^ ne phase of the 
input state by —<f>io/2 if x\ = or by +0io/2 if x\ = 1. The 
gate — [(j>oi]^ ac ^ s correspondingly on X2- The first CNOT 
gate generates the state \xi, xi(&X2) without further changing 
the phase. The line X2 contains now the result of x\ © X2 
(where © is the addition modulo 2). Subsequent application 

of — [(^J — on this line results in a conditional shift of the 

phase by — 0n/2 if x\ © X2 — or by 0n/2 if x\ © X2 = 1. 
The second CNOT restores the state |xi,X2) so that the input 
state is left unchanged apart from a phase shift that depends 
on the value of xi, X2 and x\ © X2- 



By introducing the inner product modulo two of the 
binary strings x = (xi,..., x N ) and y = (j/i, . . . , J/jv) 
(xj,yj £ {0, 1}) as x • y := X\ ■ yi © . . . © x N ■ y N we see 
that the phase shift # x according to Eq. (0 is equal to 
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(3) 



That is, the phase shifts 9 = (#oo, ■ ■ • , On) are related to 
the one-qubit rotation angles <ft = (4>ao, ■ ■ ■ , <f>ii) by an 
(unnormalized) two-qubit Hadamard transformation Ji 
(the matrix of this transformation has the entries H x , y = 
(— l) x y ). Each term in the sum of Eq. J2J represents a 
conditional phase shift by an angle <p y where the binary 
representation of the index y indicates the digits Xj that 
are part of the corresponding XOR control condition. 

Thus we can summarize that the two-qubit network 
acts on states of the computational basis by shifting their 
phase. The total phase shift is obtained by subsequently 
applying all possible conditional phase shifts which can 
be derived by combining the digits Xj of the input state. 
The conditional phase shifts are realized by applying a 
single-qubit rotation to the qubit that contains the result 
of the corresponding control condition. 

It is evident that this scheme can be generalized to an 
arbitrary number N of qubits: 

Provided that we are able to construct an 7V-qubit net- 
work which, on application to any state |x) of the N- 
qubit computational basis, 

[Nl] allows to generate — following classical logics — all 
possible control conditions x • y from the digits of 
x (where y is a binary string of length N) , 

[N2] applies z rotations to each condition x • y, 

we can implement any generalized phase shift operator 



Ug : |x) i * 



-i0 x 



|x) by using one and the same network, 



merely by adjusting the angles </> = (</>o, . . . , 4>2 N -i) of 
one-qubit z rotations. By virtue of TL 2 = 2 W 1, the angles 
(j) are related to the desired phases 9 = (#o, . . . , 9 2 n_ 1 ) 
simply by an iV-bit Hadamard transformation 
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The network can be viewed as a "black box" with 2^ — 1 
"knobs" whose settings determine the gate U 3 realized by 
the black box. By choosing (j), any phase configuration 9 
can be programmed into the network. 

The remaining task is to show that it is indeed pos- 
sible to construct such programmable networks for an 
arbitrary number of qubits. Since for many physi- 
cal implementations coupling between arbitrary qubits 
is not available or falls off with increasing distance or 
qubit number, we will assume only nearest-neighbor cou- 
pling between the qubits |25| . Further, we will show 
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FIG. 3: The three-bit network for a system with nearest-neighbor XY coupling. The dashed lines indicate the parts of the 
two-bit network which have been replaced according to rule (ii). 



that efficient network design is possible even if CNOT is 
not the natural two-qubit gate for a given implementa- 
tion HHHi. 

For the recursive method it is convenient that all con- 
trolled phase shifts appear on the last (i.e., the iV-th) 
qubit. Then the network for N + 1 qubits can be ob- 
tained from the network for N bits simply by (i) adding 
the iV+lst line with a z rotation (i.e., a phase shift con- 
trolled by xn+i), and (ii) by replacing each conditional z 
rotation on the N-th qubit in the A7-bit network by two 
rotations of the iV+lst qubit, as shown in Fig. [5] By fol- 
lowing these rules, we get z rotations for all possible XOR 
conditions of the N + 1 qubits. The conditions for the 
rotations characterized by the indices y = (2/1, ... , 2/jv+i) 
appear in their natural order if read as a binary code. 

To illustrate this, we explain the method for up to three 
qubits with an XY interaction as, e. g., for Josephson 
charge qubits coupled by SQUID loops jlalilj . The con- 
struction is analogous for other two-bit operations used 
as the basic element. The XY interaction leads to a par- 
ticularly compact design of the circuit. This is because 
the natural gate for this interaction is again a classical 
gate, namely the product of CNOT and SWAP denoted by 
CNS [13 



CNS EE 




classical computer). An example is the four-qubit net- 
work in Fig. 0] 

By realizing such a network with a physical system 
(and supplementing it with one-qubit operations such 
as, e.g., the Hadamard gate), both the Deutsch-Jozsa 
algorithm and Grover's search algorithm can be imple- 
mented. The refined version of the Deutsch-Jozsa algo- 
rithm 0, H3| requires the gate Uf : |x) 1— > e -OTr -^ x )|x) 
where / is a constant or balanced Boolean function. 
Analogously, for Grover's algorithm the gate U g has to 
be implemented with a Boolean function g: <?(x*) = 1 for 
one particular x* and g(x ^ x*) = 0. For the general- 
ized Grover algorithm |23| . several items may be marked. 
Hence the required phase shifts 8 X are given by 7r/(x) or 
tt<7(x), respectively. The corresponding one-bit rotation 
angles 4> are readily obtained by Hadamard-transforming 
9 according to Eq. J3J . 

We mention that there are other interesting applica- 
tions for these programmable networks. For example, 
by Hadamard-transforming the N-th qubit a network 
is obtained which is capable of generating a generalized 
controlled-NOT operation on the N-th qubit with an ar- 
bitrary control condition composed from the other N — 1 
qubits. Special cases are the Toffoli gate (with an arbi- 
trary number of control bits |lj]) and the carry opera- 
tion that appears in the network for the quantum adder 
used, e.g., in Shor's factoring algorithm [24| . 

Thus, the new type of network presented here offers 



The "one-qubit network" is just a single z rotation 
Xj — [(^i]^ ■ For two qubits we modify the network in 
Fig- d such that all z rotations appear on the second line 
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X 2 -ftJr^ ^.oJ^' -L^n]^-' 

From here, the derivation of the three-qubit network is 
straightforward. The result is shown in Fig. 

We emphasize that the central issue is the concept of 
these networks expressed in the requirements [Nl], [N2]. 
Since [Nl] refers only to classical logics, optimized cir- 
cuits can be found, e.g., by an exhaustive search (on a 
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FIG. 2: Recursive extension of the network. By a SWAP with 
the iV+lst qubit the phase shift <j} yi ...y N — > y:L ... yN o (con- 
trolled by the original condition) appears now on the JV+lst 
qubit. A second SWAP and a subsequent CNOT from the N- 
th to the AT + 1st qubit give the original condition xORed with 
zjv+i — this condition controls the second 2 rotation. We men- 
tion the identity = 
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FIG. 4: Optimized four-bit network for qubits with nearest-neighbor XY coupling. With this network (plus the Hadamard 
gate for each qubit), both the Deutsch-Jozsa algorithm and database search can be realized, e.g., with Josephson charge 
qubits. We have assumed qubits coupled in a chain with periodic boundary conditions. A particularly interesting feature of 
this network is that nearest-neighbor coupling allows for parallel execution of operations. Note that the z rotations and the 
two-qubit operations appear in well-separated blocks. The optimization criterion here was to minimize the number of two-qubit 
blocks. 



efficient practical solutions for a wide range of compu- 
tational tasks; this is comprehensible for modest qubit 
numbers (on the order of 10) just by comparing with ex- 
isting solutions. On the other hand, the limit of large 
N requires more careful investigation. Although the par- 
allel execution of operations (as shown in Fig. [3} can 
reduce, in principle, the complexity to 0(2 /N) it is not 
clear whether further substantial improvement is pos- 
sible. There is exponential scaling also for another re- 
source. According to Eq. Q the accuracy required for 
the one-qubit z rotations will scale ~ 2~ N . It remains 
subject for future work whether this limit can be soft- 
ened, e.g., for certain subclasses of gates Ug. 
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qubits, the networks can be constucted efficiently by or- 
dering the </>y's such that the y's form a Gray sequence 
(see Ref. []j). 



